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Abstract 

A general form for the first-order representation of the continuous, second-order linear 
structural dynamics equations is introduced in order to derive a corresponding form of 
first-order continuous Kalman filtering equations. Time integration of the resulting first- 
order Kalman filtering equations is carried out via a set of linear multistep integration 
formulas. It is shown that a judicious combined selection of computational paths and the 
undetermined matrices introduced in the general form of the first-order linear structural 
systems leads to a class of second-order discrete Kalman filtering equations involving only 
symmetric, sparse N x N solution matrices. The present integration procedure thus over- 
comes the difficulty in resolving the difference between the time derivative of the estimated 
displacement vector ( Jjx) and the estimated velocity vector (x) that axe encountered when 
one attempts first to eliminate (x) in order to form an equivalent set of second-order fil- 
tering equations in terms of (-^x). A partitioned solution procedure is then employed to 
exploit matrix symmetry and sparsity of the original second-order structural systems, thus 
realizing substantial computational simplicity heretofore thought difficult to achieve. 

+ An earlier version of the present paper without numerical experiments was presented at 
the AIAA Guidance and Control Conference, Port aland, Ore., 20-22 August 1990, Paper 
No. AIAA 90-3387. 
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Introduction 

Current practice in the design, modeling and analysis of flexible large space structures 
is by and large based on the finite element method and the associated software. The 
resulting discrete equations of motion for structures, both in terms of physical coordi- 
nates and of modal coordinates, are expressed in a second-order form. As a result, the 
structural engineering co mmuni ty has been investing a considerable amount of research 
and development resources to develop computer-oriented discrete modeling tools, analysis 
methods and interface capabilities with design synthesis procedures; all of these exploiting 
the characteristics, of second-order models. 

On the other hand, modern linear control theory has its roots firmly in a first-order form 
of the governing differential equations, e.g., (Kwakemaak and Sivan, 1972). Thus, several 
investigators have addressed the issues of interfacing second-order structural systems and 
control theory based on the first-order form (Hughes and Skelton, 1980; Arnold and Laub, 
1984; Bender and Laub, 1985; Oshman, Inman and Laub, 1987; Belvin and Park, 1989,, 
1990). As a result of these studies, it has become straightforward for one to synthesize 
non-observer based control laws within the framework of a first-order control theory and 
then to recast the resulting control laws in terms of the second-order structural systems. 

Unfortunately, controllers based on a first-order observer are difficult to express in a pure 
second-order form because the first-order observer implicitly incorporates an additional 
filter equation (Belvin and Park, 1989). However a recent work ( Juang and Maghami, 1990) 
has enabled the first-order observer gain matrices to be synthesized using only second-order 
equations. To complement the second-order gain synthesis, the objective of the present 
paper is to develop a second-order based simulation procedure for first-order obsen 3rs. 
The particular class of first-order observers chosen for study are the Kalman Filter based 
state estimators as applied to second-order structural systems. The procedure permits 
simulation of first-order observers with nearly the same solution procedure used for treating 
the structural dynamics equation. Hence, the reduced size of system matrices and the 
computational techniques that are tailored to sparse second-order structural systems may 
be employed. As will be shown, the procedure hinges on discrete time integration formulas 
to effectively reduce the continuous time Kalman Filter to a set of second-order difference 
equations. 

The paper first reviews of the conventional first-order representation of the continuous 
second-order structural equations of motion. An examination of the corresponding first- 
order Kalman filtering equations indicates that, due to the difference in the derivative of 
the estimated displacement (J^x) and the estimated velocity (i), transformation of the 
first-order observer into an equivalent second-order observer requires the time derivative 
of measurement data, a process not recommended for practical implementation. 


2 


Next, a transformation via a generalized momentum is introduced to recast the structural 
equations of motion in a general first-order setting. It is shown that discrete time numerical 
integration followed by reduction of the resulting difference equations circumvents the need 
for the time derivative of measurements to solve Kalman filtering equations in a second- 
order framework. Hence, the Kalman filter equations can be solved using a second-order 
solution software package. 

Subsequently, computer implementation aspects of the present second-order observer are 
presented. Several computational paths Eire discussed in the context of discrete said con- 
tinuous time simulation. For continuous time simulation, an equation augmentation is 
introduced to exploit the symmetry and sparcity of the attendant matrices by maintain- 
ing state dependant control Eind observer terms on the right-hand-side (RHS) of the filter 
equations. In addition, the computational efficiency of the present second order observer 
as compared to the first order observer is presented. 

Continuous Formulation of Observers 
for Structural Systems 

Linear, second-order discrete structural models can be expressed as 

Mx + Dx + Kx — Bu + Gw , x(0) = x 0 , x(0) = x 0 (1) 

u = —Z\X — Z 2 x 

with the associated measurements 

z = H\x + H 2 x + v * (2) 

where M,D,K are the mass damping and stiffness matrices of size (N x N)\ x is the 
structural displacement vector, ( N x 1); u is the active control force (m x 1); B is a 
constant force distribution matrix ( N x m); z is a set of measurements (r x 1); Hi and H 2 
are the measurement distribution matrices (r x iV); Z\ and Z 2 are the control feedback gain 
matrices (m x N)] w and v Eire zero-meem, white Gaussian processes with their respective 
covariances Q and R) and the superscript dot designates time differentiation. In the present 
study, we will restrict ourselves to the case wherein Q and R axe uncorrelated with each 
other and the initial conditions x 0 and x 0 are also themselves jointly Gaussian with known 
means and covariances. 

The conventional representation of (1) in a first-order form is facilitated by 

{ X\ = x 

x 2 = x = ii (3) 

Mx 2 = Mx = Bu + Gw — Dx 2 — Kx\ 
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which, when cast in a first-order form, can be expressed as 

{ 


Eq = Fq + Bu + Gw, q = ( x x x 2 ) T 
z — Hq + v 


( 4 ) 


where 


Mo y ■ m-°* ~ d \ • 

MM M£} 


(5) 


It is well-known that the K a.lma.rt filtering equations (Kalman, 1961; Kalman and Bucy, 
1963) for (4) can be shown to be (Arnold and Laub, 1984): 


Eq = Fq +.Bu + EPH T R~ l z 


where 


z — z — Hq , P = 


u S T 
S L 


A J Xj 1 J X 1 

1 X 2 J 1 X j 


( 6 ) 


(7) 


in which JJ and L are positive definite matrices and the matrix P is determined by the 
Riccati equation (Kwakemaak and Sivan, 1972; Arnold and Laub, 1984) 


EPE t = FPE t + EPF t - EPH t R~'HPE t + GQG t 


( 8 ) 


The inherent difficulty of reducing the first-order Kalman filtering equations given by (6) 
to second order form can be appreciated if one attempts to write (6) in a form introduced 
in (3): 

a) i\ = x 

b) x 2 = x = £\ — L\z (9) 

c) Mx 2 = —Dx 2 — Kx i + Bu + ML 2 z 


\ 


where 


L x = (H X U + H 2 S) t R-\ L 2 = (H X S T + H 2 Lf R 


Td-1 


Note from (9b) that x 2 i\. In other words, the time derivative of the estimated dis- 
placement (z) is not the same as the estimated velocity (z); hence, x x and z 2 must be 
treated as two independent variables, an important observation somehow overlooked in 
Hashernipour and Laub (1988). 

Of course, although not practical, one can eliminate x 2 from (9). Assuming x x and x 2 are 
differentiable, differentiate (9b) and multiply both sides by M to obtain 


Mx i = Mx 2 + ML\'z 


( 10 ) 
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Substituting Mi 2 from (9c) and X 2 from (9b) in (10) yields 

Mx\ = —D(x\ — L\z) — Kx\ + Bu + ML 2 Z + ML\i (11) 

which, upon rearrangements, becomes 

Mxi + Dxi + Kx\ = Bu + ML 22 + ML\z + DL\z (12) 

There are two difficulties with the above second-order observer. First, the numerical 
solution of (12) involves the computation of xi when rate measurements are made. The 
accuracy of this computation is in general very susceptible to errors caused in numerical 
differentiation of xi. Second, and most important, the numerical evaluation of z that 
is required in (12) assumes that the derivative of measurement information is available 
which should be avoided in practice. We now present a computational procedure that 
circumvents the need for computing measurement derivatives and that enables one to 
construct observers based on the second-order models. 


Second-Order Transformation of 
Continuous Kalman Filtering Equations 

This section presents a transformation of the continuous time first-order Kalman filter to 
a discrete time set of second-order difference equations for digital implementation. The 
procedure avoids the need for measurement derivative information. In addition, the spar- 
sity and symmetry of the original mass, damping and stiffness matrices can be maintained. 
Prior to describing the numerical integration procedure, a transformation based on gener- 
alized momenta is presented which is later used to improve computational efficiency of the 
equation solution. 


Generalized Momenta 


Instead of the conventional transformation (3) of the second-order structural system (1) 
into a first-order form, let us consider the following generalized momenta (Jensen, 1974; 
Felippa and Park, 1978): 


{ a) x\ — x 

b ) X2 = AMii+Cxi 


(13) 


where A and C are constant matrices to be determined. Time differentiation of (13b) 
yields 

X 2 = AMx 1 + Cx 1 (14) 
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Substituting (1) via (13a) into (14), one obtains 


x 2 — A(Bu + Gw) — (AD — C)x\ — AKx\ 


(15) 


Finally, pairing of (13b) and (15) gives the following first-order form: 

AM OlJx,! f C -I] fzil 
AD-C l\ \i 2 J [AK 0j\x 2 j 

[a(Hu + Gu>)] ^ 

The assoriated Kalman filtering equation can be shown to be of the following form: 

!][{;]■ <") 


AM 

AD-C 


where 


L x =(HxU + H 2 S) t RT 1 , L 2 = (H\S T + H 2 L) T R 


-(U.qT 


and H\ and H 2 correspond to a modified form of measurements expressed as 


z = H\X H 2 x = HxX\ + H 2 x 2 


(18) 


where 

Hx=Hx- H 2 M~ l A~ x C , H 2 = H 2 M~ l A~ l 

Clearly, as in the conventional first-order form (9), Xi and x 2 in (17) axe now two inde- 
pendent variables. Specifically, the case of A = M -1 and C = 0 corresponds to (3) with 
x 2 = xi. However, as we shall see below, the Kalman filtering equations based on the 
generalized momenta (13) offer several computational advantages over (3). 


Numerical Integration 

At this juncture it is noted that in the previous section one first performs the elimination 
of i\ in order to obtain a second-order observer, then performs the numerical solution 
of the resulting second-order observer. This approach has the disadvantage of having to 
deal with the time derivative of measurement data. To avoid this, we will first integrate 
numerically the associated Kalman filtering equation (17). 
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The direct time integration formula we propose to employ is a mid-point version of the 
trapezoidal rule: 



(19) 


where the superscript n denotes the discrete time interval t n = nh , h is the time increment 
and 6 = h/2. 

Time discretization of (17) by (19a) at the n + 1/2 time step yields 



The above difference equations require the solution of matrix equations of 2 N variables, 
namely, in terms of the two variables x £ +1 ^ 2 and x” +1/f2 , each with a size of N. To 
reduce the above coupled equations of order 2JV into the corresponding ones of order N, 
we proceed in the following way by exploiting the nature of parametric matrices of A and 
C as introduced in (13). To this end, we write out (20) as two coupled difference equations 
as follows: 

AM(i" +1/J - x”) + <5(C r Xj +1/2 - x£ +1/2 ) 

= 6AML x z n + 1 ! 2 (21) 

(AD - C)(x" +1/2 - x?) + (x£ +1/2 - x J) + 6AKx” +1/2 

= 6(AD - C)l x z n+l ! 2 + 6L 2 z n+l12 + 6ABu n+1 ' 2 (22) 

Multiplying (22) by 8 and adding the resulting equation to (21) yields 

A(M + 6D + 6 2 K)x”+ 1/2 = (AM + 6(AD - C))x ? + 8x^ 

+{8AML X + 6 2 (AD - C)L X + 6 2 L 2 }z n+1/2 + 6 2 ABu n+1 ' 2 (23) 


Of several possible choices for matrices A and B , we will examine 

f a) A = I , C = D 
\ b) A = M~ l , C = 0 


( 24 ) 
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The choice of (24a) reduces (23) to: 

(M + 8D + 8 2 K)x" +1/2 = Mi” + $x£ + 8 2 Bu n+1 ' 2 

+6{MLi+SL 2 }z n+1/2 (25) 

so that once i” +1 ^ 2 is computed, Xj"*" 1 / 2 i g obtained from (22) rewritten as 

x% +1/2 = x% + 6g n - 8Kx” +1/2 (26) 

where 

g n = Bu n+1 ^+L 2 z n+1/2 (27) 

which is already computed in order to construct the right-hand side of (25). Hence, 
_Kx™ +1 / 2 is the only additional computation needed to obtain x” +1 ^ 2 . It is noted that 
neither any numerical differentiation nor matrix inversion is required in computing x£ +1 ^ 2 . 
This has been achieved through the introduction of the general transformation (13) and 
the particular choice of the parameter matrices given by (24a). 

On the other hand, if one chooses the conventional representation (24b), the solution of 
-n+i /2 - s Stained from (23) 


(M + 6D + S 2 K )x x f 1/2 = (M + 6D)x ” + 8Mx^ 

+8{(M + 8D)L 1 +8ML 2 }z n+1/2 + 6 2 Bu n + 1 ' 2 (28) 

Once x” +1 / 2 is obtained, x£ +1 ^ 2 can be computed either by 

x” +1/2 = (x" +1/2 - x”)/8 - L x z n+1 / 2 (29) 


which is not accurate due to the numerical differentiation to obtain x x + ^ , or by (22) 


x” +1/2 = xj + 8g n - 8M~ l Kx” +l/2 - 

M~ 1 D(x” +1/2 — x”) + 8M~ 1 DLiz n+1 f 2 (30) 

which involves two additional matrix-vector multiplications, when D ^ 0, as compared 
with the choice of A = I and C = D . . Thus (24a) is the preferred representation in a 
first-order form of the second-order structural dynamics equations (1) and is used in the 
remainder of this work. 
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Decoupling Of Difference Equations 

We have seen in the previous section, instead of solving the first-order Kalman filtering 
equations of 2 n variables for the structural dynamics systems (1), the solution of the im- 
plicit time-discrete observer equation (25) of n variables can potentially offer a substantial 
computational saving by exploiting the reduced size and sparsity of M , D and K. This 
assumes that f n+1 / 2 and u n+1 / 2 are available, which is not the case since at the n th time 
step 


u n+l/2 = _Z X x n + 112 - Z 2 x” +1/2 

(31) 

* 

f n+l/2 = z „+l/2 _ tf^+1/2 _ £ 2f "+l/2 

(32) 


requires both x" +1 ^ 2 and x£ +1/ ^ 2 even if x n+1 / 2 is assumed to be known from measurements 
or by solution of (1). Note in (32), the control gain matrices are transformed by 

Z l = Z l - Z 2 M~ 1 A~ 1 C , Z 2 = Z 2 M~ x A~ l 

-*r 

There are two distinct approaches to uncouple (25) and (26) as described in the following 
sections. 

Discrete Time Update 

Equations (31) and (32) can be approximated using 

r +1/2 ~ z n _ ^ f n _ g 2i n (33) 

u "+i/2 ~-Z x x\-Z 2 x* (34) 

This approximation leads to a discrete time update of the control force and state correction 
terms which is analogous to that which exists in experiments where a finite bandwidth of 
measurement updates occurs. For discrete time approximation, the step size h = < n+1 — t n 
should be chosen to match the time required to acquire, process and output a control 
update. 

Discrete time simulation is quite simple to implement as the control force and state cor- 
rections are treated with no approximation on the right-hand-side (RHS) of (25) and (26). 
Should continuous time simulation be required, a different approach is necessary. 

Continuous Time Update 

To simulate the system given in (25) and (26) in continuous time, strictly speaking, one 
must rearrange (25) and (26) so that the terms involving x" +1/ ^ 2 and x£ +1 ^ 2 are augmented 
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to the left-hand-side (LHS) of the equations. However, this augmentation into the solution 
matrix ( M +6D+6 2 K ) would destroy the computational advantages of the matrix sparcity 
and symmetry. Thus, a partitioned solution procedure has been developed for continuous 
time simulation as described in (Park and Belvin, 1991). The procedure, briefly outlined 
herein, maintains the control force and state correction on the RHS of the equations as 
follows. 

First, x? 1 ' 2 and x£ +1/ ^ 2 are predicted by 


*"+1/2 _ 

z i p 


= x 


i > 


*"+1/2 _ -r 

X 2p ~ Z 2 


(35) 


However, instead of direct substitution of the above predicted quantity to obtain tt£ +1 ^ 2 
and 2 p +1 ^ 2 based on (31) and (32), equation augmentations are introduced to improve the 
accuracy of Up"*’ 1 ^ 2 and 2 p +1 ^ 2 . Of several augmentation procedures that axe applicable 
to construct discrete filters for the computations of u n + 1 / 2 and z"+i/ 2 , we substitute (26) 
into (31) and (32) to obtain 


u n+i/2 _ -z x x*+ 1/2 - Z 2 (x” - 8Kx" +1/2 + 
6Bu n+1 ' 2 +6L 2 z n+1 ' 2 ) 
r+l/2 »*«+!/* -,ffi*?+l/2 — 

tf 2 (xj - 8Kx* +1/ 2 + 8Bu n + »/? + 6L 2 z n+1/2 ) 


(36) 


Rearranging the above coupled equations, one obtains 


'( I + 6Z 2 B ) 6Z 2 L 2 


6H 2 B ( I + SH 2 L 2 ) ] { *" +1 / 2 / 

-Z 2 x? - (Zi - 8Z 2 K)x" +1/2 


u n +l/2 ^ _ 




+1 /2 - H 2 x 2 - (Hi - 5H 2 AT)x” +1/2 


} 


(37) 


which corresponds to a first order filter to reduce the errors in computing x 2 = Mi + Dx. 
A second-order discrete filter for computing u and z can be obtained by differentiating u 
and z to obtain 

li = —ZiXi — Z 2 x 2 
z = i — Hiii — H 2 x 2 


{s 


(38) 


and then substituting xi and x 2 from (17). Subsequently, (19) is applied to integrate the 
equations for u and z which yields 

\I+6Z 2 B + 8 2 Z 1 M- 1 B 6(Z 2 L 2 + ZiLi-¥8ZiM- 1 L 2 ) 1 / u n+1/2 1 __ 

[ 8(H 2 B + 8HiM~ x B) I + 8Hi(Li + 8M~ l L 2 ) + 8H 2 L 2 J \ z" +1 / 2 J “ 

/«*!_«/ Zi M-'(*; - SKx" +,/2 - Dx"* 1 ' 1 ) + ZiKx" + ' /2 I f 0 ) 

tw 1 HiM- x (xZ -8Kx" +ll2 Dx” +1/2 ) + H 2 Kx” +1/2 J + \ z n+1/2 - z” J } 
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The net effects of this augmentation are to filter out the errors committed in estimating 
both x\ and x 2 . Solution of (39) for ti n+1 / 2 and f n+1 / 2 permits (25) and (26) to be solved 
in continuous time for x" +1 ^ 2 and x”" 1 " 1 ^ 2 . Subsequently, (29b) is used for x" +1 and x£ +1 . 

The preceding augmentation (39) leads to an accurate estimate of the control force and 
observer error correction at the (n+1/2) time step. Although (39) involves the solution of 
an additional algebraic equation, the equation size is relatively small ( size = number of 
actuators (m) plus the number of measurements (r) ). Thus, (39) is an efficient method 
for continuous time simulation of the Kalman filter equations provided the size of (39) is 
significantly lower than the first order form of (4). The next section discusses the relative 
efficiency of the present method and the conventional first order solution. More details on 
the equation augmentation procedure (39) may be found in Park and Belvin (1991). 

Finally, it is noted that by following a similar time discretization procedure adopted for 
computing x" +1 ^ 2 and x""*” 1 ^ 2 , the structural dyn ami cs equation (1) can be solved by 

/ (M + 8D + 8 2 K)x^ l/2 = Mi" + 8x% + 8 2 Bu n+1 ' 2 . . . 

\ X 2 +1/2 = x£ + 8Bu n + x ! 2 - 8Kx"+ 1/2 1 * 

Thus, numerical solutions of the structural dynamics equation (1) and the observer equa- 
tion (20) can be carried out within the second-order solution context, thus realizing sub- 
stantial computational simplicity compared with the solution of first-order systems of equa- 
tions (4) and the corresponding first-order observer equations (6). 

It is emphasized that the solutions of both the structural displacement x and the re- 
constructed displacement x employ the same solution matrix, (M + 8D + 8 2 K). The 
computational stability of the present procedure can be examined as investigated in Park 
(1980) and Park and Felippa (1983, 1984). The result, when applied to the present case, 
can be stated as 

£ 2 Amax < 1 (41) 

where A m ax is the maximum eigenvalue of 

(\ 2 I + \Z 2 B + Z 1 M- 1 B)y = 0 (42) 

Experience has shown that |A max | is several orders of magnitude smaller than /i m ax of 

the structural dynamics eigenvalue problem: 

fiMy = Ky (43) 

Considering that a typical explicit algorithm has its stability limit Umax ■ h < 2, the 
maximum step size allowed by (42) is in fact several orders of magnitude larger than 
allowed by any explicit algorithm. 
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Computational Efficiency 

Solution of the Kalman filtering equations in second-order form is prompted by the po- 
tential gain in computational efficiency due to the beneficial nature of matrix sparcity and 
symmetry in the solution matrix of the second-order observer equations. There is an over- 
head to be paid for the present second-order procedure, that is, the additional computations 
introduced to minimize the control force and observer error terms on the right-hand-side of 
the resulting discrete equations. The following paragraphs show the second-order solution 
is most advantageous for observer models with sparse coefficient matrices M, D and K. 

Solution of the first order Kalman filter equation (6) or the second-order form (25-26, 39) 
may be performed using a time discretization as given by (19). For linear time invariant 
(LTI) systems, the solution matrix is decomposed once and subsequently upper and lower 
triangular system solutions are performed to compute the observer state at each time step. 
Thus, the computations required at each time step result from calculation of the RHS 
and subsequent triangular system solutions. For the results that follow, the number of 
floating point operations per second (flops) are estimated for LTI systems of order 0(N). 
In addition, it is assumed that the mass, damping and stiffness matrices ( M,D and K) 
are symmetric and banded with bandwidth aN, where 0>< a < (0.5 — 277). 

The first-order Kalman filter equation (6) requires (AN 2 + 2 Nr + O(N)) flops at each 
time step. The discrete time second-order Kalman filter solution (25-26, 33-34) require 
(8a 2 IV 2 -f 2alV 2 +3Nm+ANr + 0(N)) flops and the continuous time second-order Kalman 
filter (25-26, 39) require (8a N 2 + 2 aJV 2 + 5 Nm + 6 Nr + (r + m) 2 + O(N)) flops at each 
time step. To examine the relative efficiency of the first-order and second-order forms, 
several cases are presented as follows. 

First, a worst case condition is examined whereby M, D and K are fully populated (a = 
0.5 — 277) and r = m = N. For this condition, the number of flops are: 

First Order 61V 2 4- O(N) 

'< Second Order Discrete 10 N 2 + O(N) 

k Second Order Continuous 181V 2 4- O(N) 

Thus, for non-sparse systems with large numbers of sensors and actuators relative to the 
system order, the first order Kalman filter is 300 percent more efficient than the second- 
order continuous Kalman filter solution presented herein. 

For structural systems, M, and K are almost always banded. In addition, the number 
of sensors and actuators is usually small compared to the system order N. Hence, the 
value of a for which the second-order form becomes more computationally attractive than 
the first order form must be determined. If the assumption is made that the number of 
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actuators (m) and the number of measurements (r) is proportional to the bandwidth ( 
r = m = alV), the value of a which renders the second-order solution more efficient is 
readily obtained. For the 'second-order discrete Kalman filter, when a < 0.394 the second- 
order form is more efficient. Similarly, the second-order continuous K alman filter form is 
more efficient when a < 0.279. Since a obtains values approaching 0 when a modal based 
structural representation is used with few sensors and actuators, the second-order form 
can be substantially more efficient than the classical first-order form. A more detailed 
discussion can be found in Belvin (1989). 


Implementation and Numerical Evaluations 

The second-order discrete K alman filtering equation derived in (25) and (26) have been 
implemented along with the stabilized form of the controller u and the filtered measure- 
ments z in such a way the observer computational module can be interfaced with the 
partitioned control-structure interaction simulation package developed previously (Belvin, 
1989; Belvin Park, 1991; Alvin and Park, 1991). Table 1 contrasts the present CSI simula- 
tion procedure to conventional procedures. It is emphasized that the solution procedure of 
the present second-order discrete Kalman filtering equations (25) and (26) follows exactly 
the same steps as required in the solution of symmetric, sparse structural systems (or the 
plant dynamics in the jargon of control). It is this attribute that makes the present discrete 
observer attractive from the simulation viewpoint. 

The first example is a truss beam shown in Fig. 1, consisting of 8 bays with nodes 1 and 
2 fixed for cantilevered motions. The locations of actuator and sensor applications as well 
as their directions are given in Table 2. Figures 2, 3 and 4 are the vertical displacement 
histories at node 9 for open-loop, direct output feedback, and dynamically compensated 
feedback cases, respectively. Note the effectiveness of the dynamically compensated feed- 
back case by the present second-order discrete Kalman filtering equations as compared 
with the direct output feedback cases. Figure 5 illustrates a testbed evolutionary model 
of an Earth-pointing satellite. Eighteen actuators and 18 sensors are applied to the sys- 
tem for vibration control and their locations are provided in Tables 3 and 4. Figures 6, 
7, and 8 are a representative of the responses for open-loop, direct output feedback, and 
dynamically compensated cases, respectively. Note that u, response by the dynamically 
compensated case does drift away initially even though the settling time is about the same 
as that by the direct output feedback case. However, the sensor output axe assumed to 
be noise-free in these two numerical experiemnts. Although the objective of the present 
paper is to establish the computational effectiveness of the second-order discrete Kalman 
filtering equations, we conjecture that for noise-contaminated sensor output for which one 
would apply dynamic compensated strategies, the relative control performance may turn 
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out to be the opposite. Further simulations with the present procedure should shed light 
on the performance of dynamically compensated feedback systems for large-scale systems 
as they axe computationally more feasible than heretofore possible. 

Table 5 illustrates the computational overhead associated with the direct output feedback 
vs. the use of a dynamic compensation scheme by the output present Kalman filtering 
equations. In the numerical experiments reported herein, we have relied on Matlab software 
package (Wolfram, 1988) for the synthesis of both the control law gains and the discrete 
Kalman filter gain matrices. It is seen that the use of the present second-order discrete 
Kalman filtering equations for constructing dynamically compensated control laws adds 
computational overhead, only an equivalent of open-loop transient analysis of symmetric 
sparse systems of order N instead of 2 N x 2 N dense systems. 


Summary 

The present paper has addressed the advantageous features of employing the same direct 
time integration algorithm for solving the structural dynamics equations also to integrate 
the associated continuous Kalman filtering equations. The time discretization of the re- 
sulting Kalman filtering equations is further facilitated by employing a canonical first-order 
form via a generalized momenta. When used in conjunction with the previously developed 
stabilized form of control laws (Park and Belvin, 1991), the present procedure offers a sub- 
stantial computational advantage over the solution methods based on a first-order form 
when computing with large and sparse observer models. 

Computational stability of the present solution method for the observer equation has been 
assessed based on the stability analysis result of partitioned solution procedures (Park, 
1980). To obtain a sharper estimate of the stable step size, a more rigorous computational 
stability analysis is being carried out and will be reported in the future. 
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Structure: 


o) 


where 


and 


Sensor Output: b ) 

Estimator: c) 

Control Force: d) 

Estimation Error: e) 


Mq + Dq + Kq = f + Bu + Gw 

q(0) = qo, q(0) = qo 
z = Hx + v 

Mq + Dq + Kq = f + Bu + ML27 
q(o) = 0 , q(0) = 0 

u + F2M -1 Bu = F2(M -1 p + L27) + Fiq 
7 + H v L 2l = z - - Bu) - H d 4 


Table la Partitioned Control-Structure Interaction Equations 


Structure: 




Sensor Output: b) 

Estimator: c) 


Control Force: d ) 

l Estimation Error: e) 


x = Ax + Ef + Bu + Gw 
q(0) = q 0 j q(0) - q 0 
z = Hx + v 

i = Ax + Ef + Bu + L7 
x(0) = 0 
u = -Fx 

7 = z - (H<iq + H„q) 


H = [H<i H„], L 


-{?}■ ■ -m 
■si- ■ 


o' 

~ M _1 _ 

, F = [Fi F 2 ] 
Table lb Conventional Control-Structure Interaction Equatioons 


A = 


r 0 

1 




0 

1 

s 

1 

-M -1 D 


B = 

M_1 B 
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TABLE 2a: 

Actuator Placement for Truss Example Problem 


Actuator 

Node 

Component 

1 

2 

y 

2 

18 

y 

3 

9 

y 

. 4 

9 

X 


TABLE 2b: 

Sensor Placement for Truss Example Problem 


Sensor 

Type 

Node 

Component 

1 

Rate 

2 

y 

2 

Rate 

18 

y 

3 

Rate 

9 

y 

4 

Rate 

9 

X 

5 

Position 

9 

y 

6 

Position 

9 

X 


18 









TABLE 3: 

Actuator Placement for EPS Example Problem 


Actuator 

Node 

Component 

1 

97 

X 

2 

97 

z 

3 

96 

X 

4 

96 

z 

5 

65 

y 

6 

68 

y 

' 7 

59 

y 

8 

62 

y 

9 

45 

y . 

10 

45 

z 

11 

70 

y 

12 

70 

z 

13 

95 

X 

14 

95 

y 

15 

95 

z 

16 

95 

<f>T 

17 

95 

<t>y 

18 

95 
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TABLE 4: 

Sensor Placement for EPS Example Problem 


Sensor 

Type 

Node 

Component 

1 

Rate 

97 

X 

2 

Rate 

97 

z 

3 

Rate 

96 

X 

4 

Rate 

96 

z 

5 

Rate 

65 

y 

6 

Rate 

68 

y 

7 

Rate 

59 

y 

8 

Rate 

62 

y 

9 

Rate 

45 

y 

10 

Rate 

45 

z 

11 

Rate 

70 

y * 

12 

Rate 

70 

z 

13 

Position 

95 

X 

14 

Position 

95 

y 

15 

Position 

95 

z 

16 

Position 

95 

<t>x 

17 

Position 

95 

4>y 

18 

Position 

95 

<t>Z 
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TABLE 5: 

CPU Results for ACSIS Sequential and Parallel Versions 


Model 

Problem 

Type 

Sequential 

Parallel 

54 DOF 

Transient 

4.5 


Truss 

FSFB 

9.4 



K. Filter 

13.0 

10.7 

582 DOF 

Transient 

98.6 

100.3 

EPS7 

FSFB 

190.2 

294.5 


K. Filter 

284.2 

321.5 
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Truss Model: Open Loop Transient Response 

6.0 x 10“ 4 
4.8 x 10 -4 
3.6 x 10~ 4 
2.4 x 10~ 4 
1.2 X 10~ 4 


Time, sec 

Node 9, uy 

Figure 2: Truss Transient Response 
Truss Model: Full State Feedback Response 

6.0 x 10~ 4 
4.8 x 10“ 4 
3.6 x 10~ 4 
2.4 x 10~ 4 
1.2 X 10~ 4 
0.000 

0.000 0.200 0.400 0.600 0.800 1. 

Time, sec 




Node 9, uy 

Figure 3: Truss FSFB Response 




Deflection 


Truss Model: Controlled Response with Kalman Filter 
6.0 x 10“ 4 

4.8 x 10 -4 

3.6 x 10“ 4 

2.4 x 10~ 4 

1.2 x 10“ 4 

0.000 

0.000 0.200 0.400 0.600 0.800 1.000 

Time, sec 

Node 9, uy 

Figure 4: Truss Response with Filter 
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EPS7 Model: Open Loop Transient Response 

2.0 x 10” 4 

) 

1.4 x 10“ 4 

8.0 x 10 -5 

2.0 x 10 -5 

4.0 x 10“ 5 

1.0 x 10 -4 

0.000 2.000 4.000 6.000 8.000 10 

Time, sec 

— Node 45, ux Node 45, uz 

Node 45, uy 

Figure 6: EPS Transient Response 
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EPS7 Model: Full State Feedback Response 

2.0 x 10 -4 
1.4 x 10~ 4 

8.0 x 10“ 5 

2.0 x 10~ 5 
-4.0 x 10“ 5 
-1.0 x 10“ 4 

0.000 2.000 4.000 6.000 8.000 10.000 

Time, sec 

Node 45, ux Node 45, uz 

Node 45, uy 

Figure 7: EPS FSFB Response 
Model: Controlled Response w/Kalman Filter 

2.0 x 10 -4 
1.4 x 10“ 4 

8.0 x 10~ 5 

2.0 x 10“ 5 
-4.0 x 10“ 5 
-1.0 X 10 -4 

0.000 2.000 4.000 6.000 8.000 10.000 

Time, sec 

Node 45, ux Node 45, uz 

Node 45, uy 

Figure 8: EPS Response with Filter 
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